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Abstract 

We extend and develop a method for perturbative calculations of anomalous dimensions 
and mixing matrices of leading twist conformal primary operators in conformal field theories. 
Such operators lie on the unitarity bound and hence are conserved (irreducible) in the free 
theory. The technique relies on the known pattern of breaking of the irreducibility conditions 
in the interacting theory. We relate the divergence of the conformal operators via the field 
equations to their descendants involving an extra field and accompanied by an extra power of 
the coupling constant. The ratio of the two-point functions of descendants and of their primaries 
determines the anomalous dimension, allowing us to gain an order of perturbation theory. We 
demonstrate the efficiency of the formalism on the lowest-order analysis of anomalous dimensions 
and mixing matrices which is required for two-loop calculations of the former. We compare these 
results to another method based on anomalous conformal Ward identities and constraints from 
the conformal algebra. It also permits to gain a perturbative order in computations of mixing 
matrices. We show the complete equivalence of both approaches. 
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1 Introduction 



The framework of the Wilson-Kadanoff operator product expansion pQ for correlation functions 
of field operators in quantum field theory, which determine physical observables, enormously 
facilitates the analysis of their short- (or light-cone-) distance structure. Since its discovery it 
has found a large range of applications stretching from phase transition phenomena in condensed 
matter physics to scattering amplitudes in four- dimensional gauge theories. The main advantage 
of the formalism is that it allows one to evaluate a product of field operators near the short or 
light-cone distance singularity in terms of composite operators. In an interacting field theory, the 
latter mix among each other under renormalization and acquire nontrivial anomalous dimensions. 
Their evaluation at higher orders in the coupling constant is one of the goals of perturbative field 
theory. 

The recent surge of interest in anomalous dimensions of Wilson operators in the maximally 
supersymmetric Yang-Mills theory (A/" = 4 SYM) was inspired by the gauge/string correspon- 
dence which identifies them with the energies of excitations in the dual description in terms of 
a string theory on an AdSsxS 5 background [2j. The strong/weak nature of this duality makes 
its straightforward tests difficult, since the perturbative coupling expansion windows in both 
theories do not overlap. Due to complications in the quantization of string theories on warped 
backgrounds, such direct tests would require an exact evaluation of anomalous dimensions in 
gauge theory. In practical terms, what one can realistically do is to perform multiloop compu- 
tations, then use results as initial data for the conjectured integrating of M = 4 SYM. In this 
way one may be able to determine the strong-coupling asymptotics of the anomalous dimensions 
and compare it to the string theory predictions. 

The goal of the present study is to develop a formalism for efficient multiloop calculations of 
anomalous dimensions of a certain class of Wilson operators in Af = 4 SYM, namely operators 
of leading twist. An important feature of this model is that it stays anomaly- free even quantum 
mechanically and thus preserves all classical symmetries of the Lagrangian. The central object 
of our consideration is a correlation function of two conformal primary operators of leading 
twist. The conformal symmetry of the model completely determines (up to normalization) their 
functional dependence on the space-time interval between these points with the exponent given 
by the scaling dimension of the composite operators (given by the sum of their canonical and 
anomalous dimensions). A key observation for our formalism is that twist-two conformal primary 
operators lie near the conformal unitarity bound and hence are conserved in the non-interacting 
theory. However, they acquire a non- vanishing divergence in the interacting theory. It is obtained 
by applying the field equations of motion and thus is proportional to the coupling constant. 
In conformal terms, this divergence defines a particular conformal descendant of the primary 
operator. The idea of the method is to compute the ratio of the two-point correlation functions 
of descendants and of their primaries. This ratio is proportional to the anomalous dimension 
and involves an overall factor of two powers of the coupling constant. Thus, in order to evaluate 
the anomalous dimension at order n in perturbation theory, it is sufficient to compute the two 
correlators at order n — 1 and thus gain an order of perturbation theory. 

The method we describe here has been first used in [3] to calculate the one-loop anomalous 
dimensions of some twist-two operators of low spin in M = 4 SYM and then generalized to 
arbitrary spin in [3]. In this paper we present a simpler version of [4J which does not require 

lr The literature on the subject is vast, we therefore refer the reader to comprehensive proceedings of a recent 
workshop for details, [http : //www-spht . cea . f r/ Meet ings/Rencitz2007/ agend a . php| 
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super symmetry. A similar method was applied to obtain the two- [5] and three-loop [6] anomalous 
dimension of the Konishi operator and of the twist-three operator of the BMN series. 

It is a common knowledge that conformal symmetry is broken provided that the renormal- 
ization group function (3 of the coupling constant is non-vanishing. This is a direct consequence 
of dimensional transmutation which generates an intrinsic mass scale in the theory, modifying 
the scaling behavior of correlation functions. In perturbative calculations one has to use a reg- 
ularization procedure for ultraviolet divergences to render correlation functions finite. The only 
consistent regularization method for non-Abelian gauge theories is dimensional regularization 
or its spin-off, dimensional reduction. However, neither of them preserves all space-time sym- 
metries of the regularized theory, in particular, they violate the scaling and special conformal 
boosts symmetries. This has profound consequences for the form of the correlation functions even 
when the regulator is eliminated. Namely, the non-zero (3— function in D = 4 — 2e dimensions 
is (3 £ (g) = —leg + (3(g) and it induces an anomaly in the trace of the energy- momentum tensor. 
The renormalization of the product of the renormalized energy momentum tensor and conformal 
operators generates anomalous dimensions of the latter and leads to their mixing under conformal 
boosts, as we will demonstrate in this study. Even when the four-dimensional (3— function is zero 
to all orders of perturbation theory, the conformal symmetry is violated for e ^ 0. Subtracting 
divergences and sending e — > afterwards, one generates symmetry breaking contributions to 
the dilatation operator coming from terms ~ (3 £ (g)/e. This source of the symmetry breaking is a 
peculiar feature of dimensional regularization rather than an intrinsic property of the dilatation 
operator. In other words, in gauge theories with vanishing (3 function the conformal symme- 
try breaking terms can be removed by performing a scheme transformation of the dilatation 
operator and by going over to the so-called conformal scheme. This transformation does not 
affect the eigenvalues of the dilatation operator but it does change the form of the corresponding 
eigenstates, and is required for evaluation of correlation functions to preserve their diagonality. 

Our subsequent presentation is organized as follows. In Sect. (2J we give a general introduction 
to the method of computing anomalous dimensions of leading twist operators by differentiation, 
applicable to any conformal any field theory. In Sect. [31 we begin with a one- loop calculation of 
anomalous dimensions of conformal operators in six- dimensional scalar field theory with cubic 
interaction. Then we turn to two-loop order and demonstrate that conformal symmetry is broken 
and bare conformal operators start to mix in the minimal subtraction scheme. We introduce the 
so-called conformal scheme which preserves the autonomous renormalization group equation for 
conformal operators. We demonstrate the mixing for the scalar theory making use of an appro- 
priate two-point correlation function of a conformal operator and its descendant. We explain how 
this method allows us to gain an order of perturbation theory. We then turn to M = 4 SYM and 
compute the mixing matrix. Section H] is dedicated to an alternative computation of the mixing 
matrix within a different formalism to evaluate the same quantity making use of conformal Ward 
identities and commutator constraints stemming from conformal algebra. Finally, we conclude. 

2 Anomalous dimensions by differentiation 

We are interested in renormalization properties of operators of leading twist (i.e., dimension 
minus spin) in D = 2h dimensions and will restrict ourselves to those that are build from scalar 
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fields. We are going to outline our method using a simple six-dimensional model (0 3 theory), 

Nf N f 

i=i i=i 
which is conformally invariant at the classical level. Moreover, its (3 function [7] 

m = + °(9 5 ) = k N f ( 2 - 2 ) 

can be made vanish up to order g 5 in coupling constant by choosing Nf = 4. This allows us 
to use conformal symmetry arguments up to this order of perturbation theory. For the specific 
operators that we are interested in, the treatment of M = 4 SYM requires only minor changes. 
The J\f = 4 SYM Lagrangian in Minkowski space has the fornj§ 

£at=4 = tr I - \F, V F^ + | {V^ AB ) (Vy AB ) + y 2 YM [<P AB , <P CD ][<Pab, 0cd] 

+ 2\X hA cjfV% - V2g YM X aA [(f> AB , Af ] + v^^mA^^, A|] J . (2.3) 

In the following we perform perturbative expansions with respect to the coupling g = v^c^ym- 



2.1 Preliminaries 

The composite operators under consideration are built from two field operators and are generically 
written as 

3 

Qui-m = a jk d {f11 ■ ■ ■ d^ k <p d^ k+1 ■ ■ ■ <p , (2.4) 

k=0 

where {• ■ • } stands for traceless symmetrization, and if stands for an elementary scalar field of 
canonical dimension h — 1 in D = 2h space-time dimensions. In M = 4 SYM they have leading 
twist two, in the 3 theory they have twist four. From Eq. (12. 4p we can see that for a given 
spin j there are in general j + 1 coefficients to determine (one of which is just an overall 
normalization). 

There exists a continuous series of unitary irreducible representations of the conformal group 
characterized by the conformal dimension d and the Lorentz spin j > (symmetric traceless 
tensor). Unitarity requires that 

d>2(h-l)+j. (2.5) 

When the bound is saturated, these representations become reducible and one can impose irre- 
ducibility conditions. For the operators O fil '" >lj this happens in the free field theory, where they 
possess the canonical dimension 

d 0jJ =2(h-l)+j. (2.6) 
Then the irreducibility condition implies that the Ml are conserved tensors: 

d li O^-^ = Q. (2.7) 

2 We use conventions of Ref. [8]. 
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Equations ( 12. 71) allows one to fix the coefficients djk in (12.41) up to an overall normalization. 

Before writing down the solution, however, let us introduce an efficient tool for traceless 
symmetrization of the vector indices [5]. It consists in projecting all the vector indices of a 
given tensor with a null vectoi[§ z^, z 2 = 0, thus automatically symmetrizing Lorentz indices and 
suppressing all the traces. Thus we define the projected operator 



j 



6. = omi-w ^ . . . H , = J2 a jk d k cp &~ k (p , (2.8) 

k=0 

where we introduced the notation d = z^d^. To free some indices from contractions with 
z^— vectors (e.g., in order to take a divergence), one has to differentiate with respect to the 
auxiliary vector z M , in the presence of the constraint z 2 = 0. The lowest order differential 
operator which does this is second-order and reada^ [10] 



A" = (h-l + z- d z )d% - ~z»d z ■ d z . (2.9) 

So, for instance, we may recover the operator in (12.41) from the projected one in ( 12. 8ft by 

qvi-h = n A m • • • A^Oj , (2.10) 

where N is an inessential normalization constant. The symmetrization and tracelessness of 
qw-w i n ( 12.101) are automatic due to the properties 

[A'*,A u ] = 0, A^A M = (2.11) 

of the differential operator A M . Let us now rewrite (12.81) in a bi- local form, 

d j = P j (d a ,fy<p(x a )<p(x b ), (2.12) 

where 

Pj (x, y) = J2 a jk x k y j ~ k = (x + y y Pj (^-) (2.13) 

is a homogeneous polynomial of degree j. 

Using the projection variables z^ and the differential operator A M introduced in (12. 9p . we can 
rewrite (12. 7p in the following way: 

d„A»d j = 0. (2.14) 

The advantage of (12.141) over (12. 7p now becomes evident: the condition (12.141) turns into a simple 
differential equation [11] for the polynomial Pj(x) defined in (12.131) . 



((1 -x 2 )^- (2i/ + l)ar^ + j{j + 2v)\ Pj {x) = 



(2.15) 



wher^fl v — h — 3/2. Its regular solution is the Gegenbauer polynomial Cj(x) and hence the 
operators (12121) read [El [13] 

Oj = a + d b ) J q (t—^c ) (pix a )(p{x h ) . u = h 3/2 . (2.16) 
\o a + d b 



In Minkowski space this is a light-like vector, in Euclidean space it is a complex isotropic vector. 



3 

4 Owing to [A M , z 2 ] — 2z 2 g|j7 = 0, this is compatible with the constraint z 2 = 0. 
5 For fields with a light-cone spin s the index is v = h — 3/2 + s. 
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2.2 Anselmi's trick 



Here we illustrate the main idea of our method [3] . We exploit conformal properties of quantum 
field theories supposing that the renormalization has already been done. When the regulator, 
e.g., £ in dimensional regularization, is sent to zero, one expects conformal properties to emerge. 
As it is well known, renormalization introduces scheme dependence. We define our renormal- 
ization scheme (within dimensional regularization) such that conformal primary operators have 
conformal correlation functions [13] . 

When interactions are turned on (g ^ 0), the scaling dimension of the tensors Oj becomes 
in general coupling-dependent. In the following we assume that the quantum numbers of the 
composite operators are chosen in such a manner that they form a set closed under renormaliza- 
tion. After the diagonalization of the mixing matrix, the operators have well-defined conformal 
properties. In particular, their scaling dimension can be simply written as a sum of the canonical 
d j and anomalous 7j(g 2 ) dimension, 

<lj <h,j I ',(//-')• (2.17) 

Then the unitarity bound ( 12. 5ft is no longer saturated and hence the conservation equation ( 12.141) 
turns into the non-conservation equation 

d^A^d^gkj-!. (2.18) 

The right-hand side defines a (classical) conformal descendant (a state in the infinite-dimensional 
space of the conformal UIR). In our case, where the primary (lowest weight) state is a bilinear 
operator, the descendant Kj-i is trilinear in the fields and can be calculated using the classical 
field equations following from (12. ip . 

In a conformal field theory, the form of the two-point function of conformal primary operators 
Oj of conformal dimension dj and spin j is fixed to be^| 

(djfri) O k (x 2 )} = 8 jk Cj(g) P (x 2 12 y d i (2.19) 

where 

I = Ifj,u z l z 2 i I 11 = T]^ — 2 g i x \2 = X l ~ X 2 i (2.20) 



X 
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and Cj(g) are normalization constants. In (I2.19P we have used two independent projection 
variables, z% and z 2 , for the operators at points x^ and x 2 , respectively. 

The idea now is to take the divergence at both points of (I2.19p and to replace the result by 
the descendant (J23EJI0 

d^A 1 ^A 2u (6 J (x 1 ) 6j(x 2 )) = g 2 {kj^{x x ) K 3 -i(x 2 )} . (2.21) 

Using the expressions (I2.19p . (12.201) . it is straightforward to carry out the differentiation in (I2.2ip . 
When it is done, we equate the auxiliary vectors z% = z 2 = z^, for simplicity. Then we evaluate 
the ratio to be 

^^ (Kj ^Kj t (x a )) = _ 1 . (g 2 )jij + h _ 2) [ (j + h _ 1){j + 2h _ 3) (222) 
(Oj{x 1 )Oj(x 2 )) L 

+jj(g 2 )(f + hj-2j + h-l) . 



6 This equation can be found by requiring covariance under the action of the conformal group, see, e.g., 
7 Equation (|2.2ip is to be understood for non-coincident points, x\ ^ X2, otherwise one would need to consider 
possible contact terms. These terms are irrelevant for our discussion, so we can simply ignore them. 
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Figure 1: Feynman diagrams for the evaluation of anomalous dimensions (a and b) and mixing 
matrix (c) in 3 theory. The symbol (g> stands for the conformal operator. 

In a conformal field theory, the ratio (I2.22p is considered as an exact (non-perturbative) 
expression. In practice we can calculate its left-hand side only perturbatively. It is important to 
realize the appearance of the factor g 2 in the left-hand side of (12.221) . If we want to determine 

^ = c&^ + c^ )+ -- (2 - 23) 

up to, say, order (g 2 ) n , we only need to evaluate the correlators on the left-hand side of (12.221) 
up to order {g 2 ) n ~ 1 - In other words, this method allows us to gain one perturbative order in the 
calculation of 7j(g 2 ). This simple observation is the main point in our approach. 

3 Evaluation of anomalous dimensions 

In this section, we use the general formula (12.221) to perform actual loop calculations. We start 
with a simple scalar field-theory example and then turn to maximally supersymmetric gauge 
theory, which is the focus of our study. 

3.1 Anomalous dimensions in 3 theory in six dimensions 

Let us consider conformal operators of leading twist, 

6 3 = (4 + d b y cf ( ) M^ 2 (x b ) , (3.i) 

\Oa + O h J 

which are 'flavor' non-singlets (to avoid additional mixing) and are thus closed under renormal- 
izatiorH The descendants (12 . 181) are straightforwardly evaluated and are given by 

k s - x = (4 + 4 + 4r Vi ( f + t b ~t ~! a + ° b ~° c ) x(xa)M* b )<P 2 (x c ) , (3.2) 

\d a + d b + d c d a + d b + d c I 



8 The case j = 1 corresponds to the one of the conserved U(Nf) currents, M = <9 M 0i</> 2 — 0i9 M 2 , for which 
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where 

kj-iix, y) = 3 [2Cf_\(x) - 5(1 - x)C] /2 2 (x) - 2Cf_\(y) - 5(1 + y)C]L\(y)\ (3.3) 

is now a polynomial of two variables rather than of one. 

We now utilize the ratio (12.22ft of two-point functions of primaries and descendants to evaluate 
the anomalous dimension. Since we gained a factor of g 2 , it is clear that at lowest order we have 

w h ^ %H#r^ = + h - 2)( ^' + h ~ + 2h - 3) • (3 - 4) 
(0^x1)0^x2))^ 

where the left-hand side can be evaluated by a tree-level (order g°) calculation. Depending on 
whether we want to use this formula in the 3 model or in TV = 4 SYM, we can specialize it 
to either six (h = 3) or four (h = 2) dimensions to determine the corresponding leading order 
anomalous dimensions. 

A comment on the evaluation of the two point functions that arise in (13.41) is in order here. Al- 
though it is a tree level calculation, the spin dependence presents a technical difficulty Namely, 
the composite operators Oj, Kj-i contain polynomials in the projected derivatives (see, e.g., 
( 13.31) ). and expanding them one gets a result for their two-point functions in terms of multi- 
ple sums. Therefore, it is more convenient to use the Schwinger representation for the scalar 
Euclidean propagators, 

(*,(*,)#*(*»)> = 11 ^ = H [ ^a'-V"*" , (3-5) 

so that all projected derivatives acting on propagators are essentially replaced by the a— parame- 
ters. As a consequence, instead of infinite sums one obtains integrals over the polynomials that 
appear in the definition of the conformal operators, see, e.g., (12.161) . (I3.3p . which are then evalu- 
ated using standard properties of the Gegenbauer polynomials, in particular their orthogonality 
relation ^ ^ 

dx [x(l - x)p 1/2 CJ(x)C£(z) = ^fc pj^^yj^ • (3.6) 

Let us now use (13.41) (with h — 3) in order to calculate the one-loop anomalous dimensions of 
the twist-four operators (12.161) . We need to evaluate the two correlators at order g°, see Fig. [1] 
(a). The first one fixes the normalization: 

(—2r )i +k r°° r°° 
(0^)0^X2)) = 1 ^ daaj^ df3(3(a + f3Y +k 

, C f(^|) C f(^|)e X p[- (a + /3 W 2 ]. (3.7) 

The integrals with respect to Schwinger parameters are straightforwardly evaluated using (13.51) 
and yield the following result: 

02J-2 l A \2j 

(6^)6^x2)) = 5 3k (j + l)(j + 2)(2j + 2)! (47r3)2 ^ 2j+ , • (3.8) 

It indeed possesses the orthogonal form (I2.19p expected on the basis of conformal symmetry. The 
correlation function for descendants, i.e., (Kj_i(xi)Kk-i(x2)) , is obtained in a similar fashion 
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upon evaluation of the Feynman diagram in Fig. [T] (b). Its expression differs from Eq. (13.81) only 
by a factor and thus the ratio (13.41) gives the well-known leading order approximation for the 
anomalous dimensions [T5J E] : 



(0) 



6 (j + 1)(j + 2) 



(3.9) 



3.2 Scheme ambiguities 

Up to now we have been able to avoid addressing renormalization issues because our method of 
calculating 7 J -°' ) involved only tree-level calculations. However, by going higher up in perturbation 
theory, we need to discuss how the composite operators (I2.16P are renormalized. As pointed out 
above, it is misleading to assume from the start that conformal symmetry is preserved and an 
ad hoc renormalization scheme will yield diagonal correlation functions. The Poincare invariance 
alone imposes rather weak constrains on the mixing of operators: an operator with a given spin 
j will mix under renormalization with total derivatives of lower spin operators, such as 

Oji = a^Oj (no summation). (3.10) 

The bare operators Oj k , defined in terms of bare field operators, will be our basis states to 
discuss operator mixing. The renormalized operator is obtained as a superposition given by the 
renormalization matrix Z, 

j 

6,=J2 z jkO kj . (3.11) 

fc=0 

The mixing pattern implied by Poincare symmetry results in a lower triangular matrix Z, i.e., 
Z*jk = for k > j. Our goal is to define a conformal scheme, that is, to choose Z such that the 
renormalized operators Oj have conformal two-point functions (cf. (I2.19P ) 

(Q.-On) 6 k (x 2 )) = 5, k C 3 {g) P {x\ 2 )- d ^ . (3.12) 

In Eq. (I3.12p it is understood that the regulator has been removed, e — > in dimensional 
regularization. 

We determine the Z matrix from a calculation in the MS scheme and then perform an addi- 
tional finite scheme transformation. The rotation matrix is governed by the form (I3.12p of the 
two-point correlation function. In the MS scheme we define the renormalized operator insertion 

j 

[6 3 ]=Y J Z, k d kj (3.13) 

A:=0 

in terms of a renormalization matrix Z. Perturbatively, it is given by the Laurent series: 

oo 2n n y\m\(n) 

n=l ^ ' m=l 

The anomalous dimension matrix is obtained directly from a scale variation applied to Eq. (I3.13p . 
Since the bare operator does not depend on the renormalization scale fi we get the standard 
relation 

ljk{g 2 ) = - Jim /i— (\nZ) jk (g 2 ) (3.15) 
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for computation of anomalous dimension matrix from the renormalization Z-matrix. The per- 
turbative expansion of the anomalous dimension matrix is analogous to that in Eq. (12.231) . Let 
us also recall that in this scheme it is entirely determined by the residue of the Z-matrix, 

7i*U 2 )=^SV), (3-16) 

while all higher order poles are fixed from the renormalizability of the composite operators. For 
instance, up to order g A we have (for a vanishing f3 function): 

z jk - d ik-l 3 , 4 jk - d jk - [pfj J , Z jk - - ljk . (3.17) 

Beyond leading order approximation, the anomalous dimension matrix in the scheme is non- 
diagonal and has a triangular form. Note that the eigenvalues of the anomalous dimension 
matrix are given by the diagonal entries and coincide with the scale dimensions of the conformal 
operators in the conformal scheme. The scheme transformation to the latter from MS is given 
by 

3 3 

% = E B Ik[Ojk] , Zjk = £ Bj^Z mk . (3.18) 

k=0 m=k 

The finite renormalization matrix B admits the perturbative series representation 

00 2n 
n=l ^ ' 

with the expansion coefficients B^} being triangular matrices. In a conformal field theory where 
the (3— function vanishes, the anomalous dimension matrices in the two schemes are simply related 
by 

lj (g 2 )5 3k =(B- 1 1 B) jk (g 2 ). (3.20) 

Hence, the B— matrix diagonalizes the anomalous dimension matrix evaluated in the MS scheme. 
In particular, to the first nontrivial order at which the mixing phenomena occurs we have: 

T&' — CtT -<K>- (3.21) 

Note that the knowledge of , and B^ k is sufficient to reconstruct the off-diagonal entries 7-^ 
up to order g 4 . 

Let us comment on the orthogonality of conformal operators in non-integer space-time di- 
mensions. According to Eq. (12.161) . the index of the Gegenbauer polynomials is shifted by —e. 
In practical calculations, we find it convenient to use conformal covariance of the bare operator 
in 2h dimensions and therefore we define an "e— deformed" basis 

Of = (d a + d b y Cf- £ f^f J o l( ,,,wYoJ . (3.22) 

V d a + d b 
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Then instead of (13.101) we have 

% = E4A ( t ( 3 - 23 ) 

fc=0 

The practical advantage of the basis states (13.221) is that the two-point functions at order g° are 
orthogonal even in the presence of the regulator e, 

{Of (x 1 )6f(x 2 ))\ g0 ^0, j^k. (3.24) 

It is easy to recover the mixing matrix in the standard basis from the e-deformed one by expanding 
the Gegenbauer polynomials with respect to their index. We define 

§-Cr P (x)\ P= o = -2j2d v jk C" k (x) , (3.25) 

" k=0 

and have for j > k the entries 

^- (u( -"'"» (HH4--r (3 ' 26) 



Then the desired relation reads 



B$ = B'V-d%-^(~,f-^). (3.27) 
3.3 Operator mixing: perturbative calculations 

As we explained in the previous section, the conformal renormalization scheme is defined by 
the requirement that the correctly renormalized conformal primary operator Oj should lead to 
the diagonal two-point correlation functions (I3.12p . This fixes the renormalization matrix Z 
or, equivalently, Z . In this section we show how to determine B ^ (which, due to (13.271) . is 
equivalent to finding B' 1 '), by a calculation at order g 1 . 

For this purpose it is sufficient to employ Eq. (I3.12p for j > k, 

(6 j {x 1 )6 k {x 2 )) = 0, ]>k. (3.28) 
Let us now expand (13.281) up to order g 2 in the coupling constant. With the help of the expansion 

®> = 6 ? + j£rb? >6 '' ~ (2^ t W + off 1 ) • < 3 - 29 > 

we obtain 

(OM) O k (x 2 ))\ g i - ^^{^(n)^))^ = 0, j > k . (3.30) 

Here we have used the orthogonality relation (I3.24p of the operators Of at tree level, which is 
why the 1/e terms from (13.291) have also disappeared. The finite part of (13.301) determines in 
principle B 1 . However, we can gain an order in g by taking the divergence A M <9 M at point x\. 

2 

■ F b'V (d ll A> t &- k d h (x 1 ) O k (x 2 ))\ g0 = g{K j _ l {x l ) O k (x 2 ))\ g i , j > k . (3.31) 



(27r) h d 
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In Eq. (13.311) the order 0(g) correlator has still to be regularized but it turns out to be finite 
(up to a contact term) when the regulator is set to zero, see Eq. (13.341) . Notice that we cannot 
gain a further order in g by taking the divergence at the second point in (13.31 j) because then the 
equation would be trivially satisfied. 

We are going to determine B'^ from Eq. (13.311) . On the right-hand side of this equation 
we need to evaluate (Kj-i(xi) Ok(x2))\ g for j > k, see Fig. [T] (c). It turns out that the same 
calculation, but for j < k and j = k provides a useful consistency check. Firstly, the triangularity 
of the matrix B requires 

(k^ 1 (x 1 )O k (x 2 ))\ g = 0, j<k. (3.32) 
Secondly, for j = k, it follows from (13. 12[) that 

^k^d^U = (o) - ( • + h _ 2) (3 _3 3) 

Of course, we could also have taken another divergence at the second point of the correlator, as in 
(13.41) . to gain a further order in g, but our goal here is to provide a check for the calculation of the 
mixing matrix. We have indeed verified that (13.331) reproduces the previously determined values 
(13.91) . (I3.4ip of in the 3 theory and in J\f = 4 SYM, respectively, and that the triangularity 
condition (I3.32j) is satisfied. 

The calculation of (Kj_i(xi) Ok(x 2 )) \ g involves two Feynman diagrams, see Fig. [Q(c), with 
just one interaction vertex in x-space. The only Feynman integral that arises, 



d6 ~ 2£ ^ ? 2 iZlw 2 ^ = ~o7~2~TT + °( g ) . ^1^2 (3.34) 



(^ 3 ) 4 - 2£ (^ 3 ) 2 - £ 3(x? 2 ) 

is finite in dimensional regularization, up to a contact term. Recall that we always keep x 12 ^ 
and so we can drop contact terms. Taking the ratio of (Kj^i(xi) Ofc(x 2 )) \ g and (Oj(xi) Ok(x 2 )) \ g o 
according to Eq. (13.311) . we find 

B'S } =^ } df. (3.35) 
Switching back to the undeformed basis using (13.271) . we have the result : 

Bj? = "%V ■ (3.36) 

which coincides with both the explicit evaluation [7] and the conformal symmetry predictions 
US], see Eq. (13T2TD . 

Historically, there were attempts to determine the form of the renormalized operator insertion 
by a shift of the canonical dimensions in the Gegenbauer polynomials [TTl [7] . We find that the 
renormalized conformal operator insertion can be written to this order as 

1 + tJ^) 0a + d b ) j cf- £+ ^ )/2 (f^j) <PiM¥(x b ) + 0(g±) . (3.37) 

The result we obtain in 3 theory, starting with a conformal invariant operator in 6 — 2e dimen- 
sions, indeed exhibits the expectation that the index v = h — 3/2 of the Gegenbauer polynomials 
will be shifted by the amount of the anomalous dimensions, however, it contains explicitly the 
regularization parameter contrary to previous proposals. In gauge field theories such a simple 
recipe does not work [18], as we will demonstrate once more in the next section. 



11 




Figure 2: Feynman diagrams for the evaluation of anomalous dimensions in maximally super- 
symmetric gauge theory 

3.4 Application of the method to M = 4 SYM 

In M = 4 SYM one builds twist-two operators by taking bi-linear combinations of the elementary 
fields (J)ab, X A , X a and F^ v . There are many different possibilities for constructing such operators 
[SJ [19], and in general one faces an additional mixing problem between scalars, fermions and 
gluons. This happens for instance for the superconformal primary operators constructed in 
[2"Ul H]. There the mixing problem at tree level was resolved by exploiting a super-conservation 
condition required by superconformal symmetry. Here we avoid this additional mixing problem 
by considering a different member of this superconformal multiplet, which is a conformal primary 
operators in the 20' of SU (4) [Hi [21] (in which we choose the highest-weight projection of the 
<S77(4) indices of the scalar fields cj) AB ): 

Oj = (d a + d b ) j Cf ( ] tr0 12 (a)0 12 (b) . (3.38) 

\Oa + O b J 

Note that j has to be even in (13.381) . and that j = is a special case, for which Oj=o is itself the 
superconformal primary operator of the protected energy-momentum supermultiplet. For j ^ 0, 
the Oj can be obtained by acting with four supercharges [S] on the superconformal primary 
operators of [20], H] . The fact that these operators belong to the same supermultiplet implies that 
their anomalous dimensions are given by the same universal formula. 

The expression of the descendants Kj is given in the Appendix, see Eq. (lA.lj) . Schematically, 
they read 

= A 3 ai d b , 4) tr ({A£(ar a ), A a4 (x c )}0 12 (x 6 ) + {X x & {x a ), A" 2 (x c )}0 12 (x 6 )) 

+ B 3 (d a ,d b ,d c ) tr [d^ 12 (x a )^ 12 (x b )}F^(x c )z u + 0(g) , (3.39) 

Let us first calculate 7 J -°' ) from (13. 4p . The calculations presented here were done in the light cone 
gauge, i.e., z ■ A = 0. However, because the correlators we calculate are gauge invariant, our 
results do not depend on the gauge choice. We have explicitly checked that a calculation in the 
standard covariant gauge leads to an identical result. For the (Kj-i(x\) Kj^i(x 2 )) correlator, 
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there are two contributing Feynman diagrams, one involving fermions in Fig. [2] (b) and the other 
one gluons in see Fig. [2] (c). Their respective contributions to 7] are 

if ) big. 2(a) = 2 , 7 f |Fig. 2(b) = 2 ^ - 1) , (3.40) 

where Sj = Y^k=i ^/k is a harmonic sum. Adding the two, we obtain the well-known result 

7] 0) = 2S 3 . (3.41) 

Now turning to the mixing matrix, in the left-hand side of (13.311) . the prefactor of can 
be easily evaluated by a tree-level calculation, 

%A llM d i - k (d h (x 1 )d k (x 2 )) !fl = --^M jk (x 12 ), (3.42) 

where 

(N 2 -1) 1 ( x 2 \ j+k ~ 1 
M jk (x 12 ) = 2(j + k + iy-j^f- & (-2^) , (3-43) 

and dj k is defined in (13.261) . In order to evaluate the right-hand side of (I3.3ip we need to compute 
the correlator (Kj^iOj)\ g (see Fig. [3]). We obtain 

(i^-i(xi) 6 k {x 2 )) = M jk (x) (I jk + J jk ) , (3.44) 

where the integrals Ij k , Jj k come from Fig. [3] (a) and the sum of Figs. [3] (b,c), respectively. They 
are defined by (for j, k even) 

-1 

4/2 



and 



dt (1 + t) aj (t)C L k ,z (t) = -2 (j > k) , (3.45) 
1 



J jk = j ds j dt{l+t)b j {s,t)C 1 k ,2 {s)+ j dt^-t J ds (1- s)bj(s,t) Cl /2 (t). (3.46) 

The expressions for aj(t) and bj(s,t) can be found in ( 1A.3I) . With the method described in the 
appendix of [22J one can evaluate Jj k , the result being 

J jk = 2 faSj-k + S^-Si^- 2Sj + l) (j > k) . (3.47) 

Combining the two contributions (13.451) . (I3.47P and using (13.311) . we obtain for the mixing matrix 

B'jk = 2d ]k fa +S^- 25 i _ fc - Si±k) . (3.48) 

Taking into account Eq. (13.271) and 7j°' ) = 2Sj we finally arrive at 

B fk = 2d){ 2 (Sj + S k + S^-S^- 2S' i _ fc ) . (3.49) 

This result constitutes one of the main points of our paper. Namely, we have been able to 
determine the first correction to the mixing matrix of a particular type of twist-two operators in 
the M = 4 SYM by just performing an order g 1 perturbative calculation. The standard quantum 
field theory approach would require going to order g 4 to obtain the same result. This clearly 
demonstrates the power of conformal invariance applied to higher-loop calculations. 
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Figure 3: Feynman diagrams needed for the evaluation of the mixing matrix in M = 4 SYM. 

4 Origin of the mixing matrix 

As we established in the preceding sections, even in a theory with a vanishing (3— function, 
conformal symmetry is broken in the MS scheme and leads to a mixing of conformal composite 
operators under scale transformations. To analyze this mechanism in greater detail, we employ 
conformal Ward identities [23j 122]. They can be simply derived form the reparameterization 
invariance of the generating functional, given as a path integral. To derive the true Ward identities 
it is crucial that the action contains a regulator, which in our case is done by changing the 
integration volume, i.e., d 4 x — > d D x and replacing the coupling g —> ^ D ~ 4 ^ 2 g. Performing 
an infinitesimal field transformation — > = + 8&(x) yields for a generic Green 

function of a (composite) operator 0($) to the Ward identities, 

(0($)6 (Offa) ■ ■ ■ $(x n ))> = -((50($)) $(xi) • • ■ $(£„,)> " (SiS($)) • • ■ 

(4.1) 

in the regularized theory. In a second step, we will perform the renormalization procedure, 
which yields anomalous terms that emerge from contact terms in the product of the symmetry 
variation of the regularized action functional and the operator insertion, i.e., 0($) (SiS (<&)). In 
particular when the conformal variation is involved, it generates the anomalous dimensions and 
the mixing matrix. However, the source of the latter symmetry breaking is a peculiar feature of 
dimensional regularization. In a third step, we employ conformal constraints, arising from the 
conformal algebra, to perform a scheme transformation that restores conformal symmetry such 
that conformal operators obey autonomous renormalization group equations. 

Let us first with the conformal variations of the dimensional regularized M = 4 SYM action 
(12.31) . Dilatation and special conformal boost variations lead to the following anomalous operator 
insertions: 

5 D S = eJ d D x!^0 Ak (x) + B -n^x)-n x - x (x)^ , (4.2) 
6%S = e j d D x2x^l^O Ak (x) + B -n^x)-n x - x (x)^ -2(D-2) j d D x 0^ B {x) .(4.3) 
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Here the relevant gauge invariant operator insertions, of type-A in terminology of Ref. [24]. are 

Al (x) = tvF^F^ix), (4.4) 

O m {x) = 9 \\<P ab {x)A CD (x)}[^ ab (x)Acd{x)], (4.5) 

O m (x) = V2gtT{hA(x)[<P AB {x),XUx)}-X aA {x)[^ AB (x),X B (x)}} . (4.6) 

For the composite operator (I3.38j) considered here, only the first two operator insertions Ai 
can potentially contribute to our leading order analyses of conformal anomalies. The type-B 
operators are BRST exact variations, and are irrelevant for the present study. Finally Q are 
equations of motion operator, e.g., 

^(*) = ^0 AB (*). (4-7) 

In the MS scheme, the scale and special conformal Ward-identities for the Green functions 
with the conformal operator Oj = Oj(x = 0) insertion, see Eq. (13.381) . can be found by plugging 
in the variations (14.21) and (14.31) into the generic Ward identity (14.11) . 

{[6 3 ]5 D X) = -((5 D [6 3 })X) - ([6 3 ](8 D iS)X) (4.8) 
([6 3 ]6 K X) = -(mdj))*) - ([6 3 }(5^S)X) . (4.9) 

Without loss of generality, we specify the field monomial of elementary fields to the two scalar 
fields of the theory X = <p AB (x\)(j) CD (22) • By definition, the left-hand side of these equations 
are finite, however, the separate terms on the right-hand side contain anomalous contributions 
or even divergencies. To obtain the renormalized Ward identities we will now renormalize the 
operator product of composite insertion and conformal variation of the action and finally remove 
the regularization by sending e to zero. 

In the case of the dilatation Ward identity (I4.8p . the variation in the first term on the right- 
hand side leads to the canonical dimension of the operator, while the second term in the right- 
hand side of Eq. (14.81) is ill-defined since it involves an operator product at coincident space-time 
points, 

[6j](5 D S) = e [ d D xO Al (x)[dj] - [ d D xn^(x)[dj] + ■■■ . (4.10) 



Hence a subtractive renormalization procedure is required, where the divergence is a \je— pole 
and its residue is nothing but the anomalous dimension of the operator. The pole will be cancelled 
by the e— term, appearing in the variation of the action (14.21) . A rigorous treatment can be done 
by means of differential operator vertex insertions [25J. As expected, it can be demonstrated in a 
straightforward manner that finally the Ward-identity (14.81) turns into the renormalization group 
equation for Green function with renormalized composite operator insertions, 



.4 + 27 * (fl2) 



(Pi]X) = -j2l jt (g 2 )(ld jt ]X) . (4.11) 



k=0 



Note that a (3— function proportional term is absent in our conformal theory, nevertheless, the 
trace anomaly in the regularized theory turns into the anomalous dimensions and, as we will now 
see, it is also responsible for the mixing of operators. 
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Figure 4: The symbol ® stands for the conformal operator and o stands for the operator 
insertions Oa ■ 

The renormalization procedure for the right-hand side of the special conformal Ward identity 
(14.91) has two peculiarities. First, the variation of the renormalized operator insertion, i.e., 8_ [Oj] 
leads to an infinite expression which has to be cancelled against a singularity arising from the 
renormalization of the operator product 



[Oj] = e J d D x2x.O Al (x)[d J ] - J d D x2x^(x)[O j ] + ■ 



(4.12) 



Second, the subtractive renormalization of the latter is not automatically fixed by the one ap- 
pearing in the dilatation variation. The renormalized Ward identity can be written as 



(i0 3l ]5 K X) = -i£ [a(0 +i c (l;g 2 )] «<Vi]*> + • • • , 



(4.13) 



k=0 



where the ellipsis contain BRST variations, the scaling dimensions in the conformal variation 
55 is now given by 1 + 7^, djk(l) = 2(1 — k)(l + k + 1)5 jk- Again in a conformal theory a /3 
proportional term is absent. The so-called special conformal anomaly 



(2tt) 2 v " Al 



2 7 ^ + 2[4 1 |,6])+0(^), 



(4.14) 



1 /2 ________ 

where bjk = 2(j — k)(j + k + l)d-' k , see Eq. ( I3.26p . is expressed by two subtractive renormalization 

constants = <5jfc(7 J -°' ) — 2 7( ^ ' ) )/2 and Z^~. While the former is nothing else as the residue 
of the renormalization matrix and is expressed in terms of the anomalous dimensions, cf. Eq. 
( I3.17p . the latter one is defined as 



d D x 2x-0 Al (x)[Oji 



d D x 2x-O Al {x)0 
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(4.15) 



fc=0 



and is straightforwardly evaluated to lowest order from the Feynman diagrams shown in Fig. HI 
The result is 



9 2 1 
{2ufe 



~{2Z [ lb-w)+0{g±) 



(4.16) 



16 



where the matrix elements of the w-matrix are given for j > k by 



w jk = 2{2k + 1) [1 + (-l) j ~ k ] (Sj + - 25 i _ fc ) . (4.17) 

Combining Eqs. ( 14. 1411 . (14.161) . and (14.171) . we find that the special conformal anomaly, 

7j fc (j; 9 2 ) = -2(3 -k)(j + k + l)j£pB§> + Otf) , (4.18) 

is expressed by means of the mixing matrix -B^jP, found earlier in Eq. (I3.49p . 

The anomalous dimensions and the so-called special conformal anomaly are not entirely in- 
dependent quantities. In the case of a conformal field theory it is easy to derive a constraint 
between them by acting with the differential operator [idjd\x on the conformal boost Ward iden- 
tity (I4.13P and with the generator of the conformal boost, — a differential operator acting on the 
arguments of elementary fields, — on the renormalization group equation (14. lip . Subtracting 
both equations yields zero on the left-hand side while the right-hand side is the desired constraint 
between the conformal anomalies, which we write as 

3 

2(7 -k)(j + k+ l)i jk {g) = hjm{g)i c mk {3\ g) - l] m (3\ g)lmk{g)] • (4.19) 

m=k 

The result can be understood as a consequence of the conformal commutator [D, KJ\ = iK^ 
applied to the Green function with conformal operator insertion. We note that the additional 
/—dependence in 7| fc (/; g) is governed by another conformal constraint that arises from the com- 
mutator relation [K-, P+] = —2i(D + M_ + ) and reads 

7 C (/ + l;g) = i c (l,g) - 2i(g) . (4.20) 

The above conformal constraints guarantee that there exist a scheme in which the covariant 
behavior of conformal operator under conformal transformations is ensured |14j . Such a scheme is 
obtained by a finite renormalization group transformation (I3.20p that diagonalizes the anomalous 
dimension matrix. Utilizing the constraints (14.191) and (I4.20p one finds the S-matrix in terms of 
the special conformal anomaly: 

B = = 1 - J7 C + Jh c J7 C ) , where (J 1 c ) ik = - , J/ k , r. (4.21) 

1 + J 7 C ' u ' ; ' v 1,3 2(j - k)(j + k + 1) K ' 

Needless to say that the mixing matrix to order g 2 follows by inserting the expression (14.181) for 
the special conformal anomaly into this equation and coincide with the result (13.491) . obtained 
in Sect. 13.41 Rotating now the conformal operator via Eq. (13.181) . the conformal boost Ward 
identities (14.131) turns into [H]: 

(OiiS^X) = -i 2(j -l)(j + l + l + 7i (/)) ([6 jt i-i]X) + .... (4.22) 

Here the ellipsis stands for Green functions with BRST-exact operator insertions Ob, which do 
not contribute in gauge invariant quantities. As we see, the conformal operators Oji transform 
covariantly under conformal boost and in particular the lowest state in the module Oj is invariant 
under conformal boost. 
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5 Conclusions 



As we demonstrated in this work, in a conformal field theory there exists a special renormaliza- 
tion scheme in which renormalized conformal operators possess diagonal two-point correlation 
functions ( 12.191) . whose fall-off with distance is determined merely by their scale dimensions. We 
proposed to use these correlators for evaluation of anomalous dimensions of conformal opera- 
tors at higher orders making use of the so-called Anselmi's trick. The main advantage of this 
formalism for multi-loop calculations of leading twist anomalous dimensions, compared to dia- 
grammatic evaluation of operator matrix elements, arises from simpler topologies of contributing 
Feynman graphs taking the form of bubble diagrams. Along this line of reasoning, the knowledge 
of the four loop diagrams, appearing in the descendants, would permit a three-loop calculation 
of anomalous dimensions. 

We finally remark that in a non-conformal theory, e.g., QCD, the covariant behavior of 
conformal operators is spoiled by a term proportional to the /3— function. Hence, the form of the 
correlation functions will change too. However, setting by hand the f3— function to zero at a given 
order of perturbation theory, the conformal prescription still applies and can be used for highly 
non-trivial predictions. We can even argue that in certain cases the trace anomaly can be entirely 
incorporated into the conformal predictions. For instance, we can introduce a scheme in which 
the conformal operators of leading twist are multiplicatively renormalizable. Then the form of 
the two-point correlation function in the full theory is entirely fixed by the renormalization group 
equation: 



where the running coupling satisfies the initial condition g(fi] g) = g. 

Another potential application of the method is the evaluation of the mixing matrix in full 
QCD to two- loop order accuracy, including the j3— function (see, e.g., Ref. [22]). This allows one 
to restore the full anomalous dimension matrix and even the non-forward evolution kernels in 
the MS scheme to three-loop level. This piece of information is required in the application of 
perturbative QCD to exclusive processes and would allow for a complete next-to-next-to-leading 
analysis in the MS scheme. 
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(A.B. and D.M.) are grateful to LAPTH (Annecy) for the warm hospitality extended to them 
at different stages of the work. J.H. acknowledges the warm hospitality extended to him by the 
Theory Group of the Dipartimento di Fisica, Universita di Roma "Tor Vergata" . 

A Appendix: Calculation of descendants 

We calculate the divergence (12.181) of the twist-two operator (13.381) on the classical level, using the 
equations of motion. The terms in the descendant K have two origins: the commutation relation 



(O^OOfcfo))^) = SjkCj (g(l/^r;g)) P (4,)"* exp 




(5.1) 
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of covariant derivatives and the use of the Klein-Gordon equation derived from the action (12.31) . 
The final result can be written as 

Kj. x = ^(4,4,4) tr ({A§ , (x a ),A tt4 (x c )}0 12 (a; 6 ) + {A^(x a ) > A d!2 (x c )}0 12 (x 6 )) 

+ B,(d a , 4, 4) tr [d^ 12 (x a ), 12 (x fe )]F^(x c )^ + 0( 5 ) , (A.l) 

where Aj and Bj are homogeneous polynomials of degree j — 1 and j — 2, respectively. Here we 
have dropped 0(g) terms, which are quadratic in <fi or contain two F^'s, since they are irrelevant 
for the considered perturbative order. In terms of new variables u = a + b + c, s — (a + c — b)/u 
and t = (a — c — b)/u they can be written as 

A(a, b, c) = 2v / 2V~ 1 a j (s) , Bj(a, b, c) = -i8u j ~ 2 bj(s, t) (A.2) 

with 

6,M) = -^k : («) + («-l)^'(«)l+7-^kl-^(«)-(l-0^Wl (A.3) 



s-t 



s-tf 



aj{s) = P J (s) + (s-l)P;(s), (A.4) 

and Pj(x) = Cj (x). Note there is no singularity in bj(s,t) for s = t, which can be seen by 
Taylor expanding around s = t. 
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